Search for: All records

Abstract. Changes in ocean temperature and salinity are expected to be an important determinant of the Greenland ice sheet's future sea level contribution. Yet, simulating the impact of these changes in continental-scale ice sheet models remains challenging due to the small scale of key physics, such as fjord circulation and plume dynamics, and poor understanding of critical processes, such as calving and submarine melting. Here we present the ocean forcing strategy for Greenland ice sheet models taking part in the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), the primary community effort to provide 21st century sea level projections for the Intergovernmental Panel on Climate Change Sixth Assessment Report. Beginning from global atmosphere–ocean general circulation models, we describe two complementary approaches to provide ocean boundary conditions for Greenland ice sheet models, termed the “retreat” and “submarine melt” implementations. The retreat implementation parameterises glacier retreat as a function of projected subglacial discharge and ocean thermal forcing, is designed to be implementable by all ice sheet models and results in retreat of around 1 and 15 km by 2100 in RCP2.6 and 8.5 scenarios, respectively. The submarine melt implementation provides estimated submarine melting only, leaving the ice sheet model to solve for the resulting calving and glacier retreat and suggests submarine melt rates will change little under RCP2.6 but will approximately triple by 2100 under RCP8.5. Both implementations have necessarily made use of simplifying assumptions and poorly constrained parameterisations and, as such, further research on submarine melting, calving and fjord–shelf exchange should remain a priority. Nevertheless, the presented framework will allow an ensemble of Greenland ice sheet models to be systematically and consistently forced by the ocean for the first time and should result in a significant improvement in projections of the Greenland ice sheet's contribution to future sea level change. 

Abstract. The ice sheet model intercomparison project for CMIP6 (ISMIP6) effort brings together the ice sheet and climate modeling communities to gain understanding of the ice sheet contribution to sea level rise. ISMIP6 conducts stand-alone ice sheet experiments that use space- and time-varying forcing derived from atmosphere–ocean coupled global climate models (AOGCMs) to reflect plausible trajectories for climate projections. The goal of this study is to recommend a subset of CMIP5 AOGCMs (three core and three targeted) to produce forcing for ISMIP6 stand-alone ice sheet simulations, based on (i) their representation of current climate near Antarctica and Greenland relative to observations and (ii) their ability to sample a diversity of projected atmosphere and ocean changes over the 21st century. The selection is performed separately for Greenland and Antarctica. Model evaluation over the historical period focuses on variables used to generate ice sheet forcing. For stage (i), we combine metrics of atmosphere and surface ocean state (annual- and seasonal-mean variables over large spatial domains) with metrics of time-mean subsurface ocean temperature biases averaged over sectors of the continental shelf. For stage (ii), we maximize the diversity of climate projections among the best-performing models. Model selection is also constrained by technical limitations, such as availability of required data from RCP2.6 and RCP8.5 projections. The selected top three CMIP5 climate models are CCSM4, MIROC-ESM-CHEM, and NorESM1-M for Antarctica and HadGEM2-ES, MIROC5, and NorESM1-M for Greenland. This model selection was designed specifically for ISMIP6 but can be adapted for other applications. 

Abstract. The effect of the North Atlantic Ocean on the Greenland Ice Sheet through submarine melting of Greenland's tidewater glacier calving fronts is thought to be a key driver of widespread glacier retreat, dynamic mass loss and sea level contribution from the ice sheet. Despite its critical importance, problems of process complexity and scale hinder efforts to represent the influence of submarine melting in ice-sheet-scale models. Here we propose parameterizing tidewater glacier terminus position as a simple linear function of submarine melting, with submarine melting in turn estimated as a function of subglacial discharge and ocean temperature. The relationship is tested, calibrated and validated using datasets of terminus position, subglacial discharge and ocean temperature covering the full ice sheet and surrounding ocean from the period 1960–2018. We demonstrate a statistically significant link between multi-decadal tidewater glacier terminus position change and submarine melting and show that the proposed parameterization has predictive power when considering a population of glaciers. An illustrative 21st century projection is considered, suggesting that tidewater glaciers in Greenland will undergo little further retreat in a low-emission RCP2.6 scenario. In contrast, a high-emission RCP8.5 scenario results in a median retreat of 4.2 km, with a quarter of tidewater glaciers experiencing retreat exceeding 10 km. Our study provides a long-term and ice-sheet-wide assessment of the sensitivity of tidewater glaciers to submarine melting and proposes a practical and empirically validated means of incorporating ocean forcing into models of the Greenland ice sheet. 

Abstract. Projection of the contribution of ice sheets to sea level change as part ofthe Coupled Model Intercomparison Project Phase 6 (CMIP6) takes the formof simulations from coupled ice sheet–climate models and stand-alone icesheet models, overseen by the Ice Sheet Model Intercomparison Project forCMIP6 (ISMIP6). This paper describes the experimental setup forprocess-based sea level change projections to be performed with stand-aloneGreenland and Antarctic ice sheet models in the context of ISMIP6. TheISMIP6 protocol relies on a suite of polar atmospheric and oceanicCMIP-based forcing for ice sheet models, in order to explore the uncertaintyin projected sea level change due to future emissions scenarios, CMIPmodels, ice sheet models, and parameterizations for ice–ocean interactions.We describe here the approach taken for defining the suite of ISMIP6stand-alone ice sheet simulations, document the experimental framework andimplementation, and present an overview of the ISMIP6 forcing to beused by participating ice sheet modeling groups. 

Abstract. Ice flow models of the Antarctic ice sheet are commonly used to simulate its future evolution inresponse to different climate scenarios and assess the mass loss that would contribute tofuture sea level rise. However, there is currently no consensus on estimates of the future massbalance of the ice sheet, primarily because of differences in the representation of physicalprocesses, forcings employed and initial states of ice sheet models. This study presentsresults from ice flow model simulations from 13 international groups focusing on the evolutionof the Antarctic ice sheet during the period 2015–2100 as part of the Ice Sheet ModelIntercomparison for CMIP6 (ISMIP6). They are forced with outputs from a subset of models from theCoupled Model Intercomparison Project Phase 5 (CMIP5), representative of the spread in climatemodel results. Simulations of the Antarctic ice sheet contribution to sea level rise in responseto increased warming during this period varies between −7.8 and 30.0 cm of sea level equivalent(SLE) under Representative ConcentrationPathway (RCP) 8.5 scenario forcing. These numbers are relative to a control experiment withconstant climate conditions and should therefore be added to the mass loss contribution underclimate conditions similar to present-day conditions over the same period. The simulated evolution of theWest Antarctic ice sheet varies widely among models, with an overall mass loss, up to 18.0 cm SLE, in response to changes in oceanic conditions. East Antarctica mass change varies between −6.1 and8.3 cm SLE in the simulations, with a significant increase in surface mass balance outweighingthe increased ice discharge under most RCP 8.5 scenario forcings. The inclusion of ice shelfcollapse, here assumed to be caused by large amounts of liquid water ponding at the surface ofice shelves, yields an additional simulated mass loss of 28 mm compared to simulations without iceshelf collapse. The largest sources of uncertainty come from the climate forcing, the ocean-induced melt rates, thecalibration of these melt rates based on oceanic conditions taken outside of ice shelf cavitiesand the ice sheet dynamic response to these oceanic changes. Results under RCP 2.6 scenario basedon two CMIP5 climate models show an additional mass loss of 0 and 3 cm of SLE on average compared tosimulations done under present-day conditions for the two CMIP5 forcings used and displaylimited mass gain in East Antarctica.